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Abstract 

We present improved measurements of the differential production rates of stable charged 
particles in hadronic Z° decays, and of charged pions, kaons and protons identified over a 
wide momentum range using the SLD Cherenkov Ring Imaging Detector. In addition to 
flavor- inclusive Z° decays, measurements are made for Z° decays into light (u, d, s), c and b 
primary flavors, selected using the upgraded Vertex Detector. Large differences between the 
flavors are observed that are qualitatively consistent with expectations based upon previously 
measured production and decay properties of heavy hadrons. These results are used to test 
the predictions of QCD in the Modified Leading Logarithm Approximation, with the ansatz 
of Local Parton-Hadron Duality, and the predictions of three models of the hadronization 
process. The light-flavor results provide improved tests of these predictions, as they do 
not include the contribution of heavy-hadron production and decay; the heavy-flavor results 
provide complementary model tests. In addition we have compared hadron and antihadron 
production in light quark (as opposed to antiquark) jets. Differences are observed at high 
momentum for all three charged hadron species, providing direct probes of leading particle 
effects, and stringent constraints on models. 
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1 Introduction 



The production of jets of hadrons from hard partons produced in high energy collisions 
is understood qualitatively in terms of three stages. Considering the process e + e~ — > qq, 
the first stage, "fragmentation", involves the radiation of gluons from the primary quark 
and antiquark, which in turn may radiate gluons or split into qq pairs until their virtuality 
approaches the hadron mass scale. This process is in principle calculable in perturbative 
QCD, and three approaches have been taken so far: i) differential cross sections have been 
calculated [1] for the production of up to 4 partons to second order in the strong coupling 
a s , and leading order calculations have been performed recently for as many as 6 partons 
(see e.g. [2]); ii) certain parton distributions have been calculated to all orders in a s in the 
Modified Leading Logarithm Approximation (MLLA) [3]; iii) "parton shower" calculations 
[4] have been implemented in Monte Carlo simulations that consist of an arbitrary number 
of q—>qg, g-^gg and g^qq branchings, with each branching probability determined from 
QCD at up to next-to-leading logarithmic order. 

In the second stage these partons transform into "primary" hadrons. This "hadroniza- 
tion" process is not understood quantitatively and there are few theoretical predictions for 
light hadrons, i.e. those that do not contain a c or b quark. One approach is to make the 
ansatz of local parton-hadron duality (LPHD) [3], that inclusive distributions of primary 
hadrons are the same, up to a normalization factor, as those for partons. Calculations using 
MLLA QCD, cut off at a virtual parton mass comparable with the mass of the hadron in 
question, have been used in combination with LPHD to predict properties of the distribution 
of £ = — \n(2p/ Eqm), where p and Eqm are the hadron momentum and the e + e~ energy, re- 
spectively, in the e + e~ center of mass frame. For a given primary hadron species: the shape 
should be described well by a Gaussian function within about one unit of the peak position; 
a wider £ range should be described with the addition of small distortion terms; and the 
peak position should depend inversely on the hadron mass and logarithmically on Eqm- At 
sufficiently high x p = 2p/E C M, perturbative QCD has also been used to calculate the E C m 
dependence of the x p distributions [5]. It is desirable to test these predictions experimentally 
and to encourage deeper theoretical understanding of the hadronization process. 

In the third stage unstable primary hadrons decay into the stable particles that traverse 
particle detectors. This stage is understood inasmuch as proper lifetimes and decay branching 
ratios have been measured for many hadron species. However, these decays complicate 
fundamental hadronization measurements because many of the stable particles are decay 
products rather than primary hadrons, and it is typically not possible to determine whether 
a given detected hadron is primary. Previous measurements at e + e~ colliders (see e.g. [6, 7]) 
indicate that decays of vector mesons, strange baryons and decuplet baryons produce roughly 
two-thirds of the stable particles; scalar and tensor mesons, and radially excited baryons have 
also been observed [7], which contribute additional secondaries. Ideally one would measure 
every possible hadron species and distinguish primary hadrons from decay products on a 
statistical basis. A body of knowledge could be assembled by reconstructing heavier and 
heavier states, and subtracting their known decay products from the measured differential 
production rates of lighter hadrons. The measurements presented here constitute a first step. 

Additional complications arise in jets initiated by heavy quarks, since the leading heavy 
hadrons carry a large fraction of the beam energy, restricting that available to other pri- 
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mary hadrons, and their decays produce a sizable fraction of the stable particles in the jet. 
Although the production and decay of some B and D hadrons, containing a b and c quark, 
respectively, have been studied inclusively [8] and the effect of heavy quark masses on the 
charged multiplicity in their jets has been observed [9, 10], there remain large uncertainties 
in, e.g., B® and heavy baryon decays, and heavy quark masses. The removal of heavy flavor 
events will therefore allow the direct study of the hadronization of light quark jets, and 
measurements of c and b jets may shed additional light on some of the above issues. 

A particularly interesting aspect of light quark hadronization is the question of what 
happens to the quark or antiquark that initiated a given jet. Recent experimental results [11, 
12, 13] have confirmed the notion that it is "contained" as a valence constituent of a particular 
hadron, and that this "leading" hadron has on average a higher momentum than the other 
hadrons of the same type in the jet. More precise measurements of such leading particle 
effects could lead to methods for identifying jets of a specific light primary flavor, with a 
number of applications in ep and hadron-hadron collisions as well as in e + e~ annihilations. 

There are several phenomenological models of jet production that combine modelling of 
all three stages of particle production and it is important to test their predictions. To model 
the parton production stage the HERWIG 5.8 [14], JETSET 7.4 [15] and UCLA 4.1 [16] event 
generators use a combination of first order matrix elements and a parton shower. To model 
the hadronization stage, the HERWIG model splits the gluons produced in the first stage into 
qq pairs, and these quarks and antiquarks are paired up locally to form colorless clusters that 
decay into the primary hadrons. The JETSET model takes a different approach, representing 
the color field between the partons by a semi-classical string; this string is broken, according 
to an iterative algorithm, into several pieces, each of which corresponds to a primary hadron. 
In the UCLA model whole events are generated according to weights derived from the phase 
space available to their final states and the relevant Clebsch-Gordan coefficients. Each of 
these models contains arbitrary parameters that control various aspects of hadronization and 
have been tuned to reproduce data from e + e~ annihilations. The JETSET model includes a 
large number of parameters that control, on average, the species of primary hadron produced 
at each string break, giving it the potential to model the observed properties of identified 
hadron species in great detail. In the HERWIG model, clusters are decayed into pairs of 
primary hadrons according to phase space, and the relative production of different hadrons 
is effectively governed by two parameters that control the distribution of cluster masses. In 
the UCLA model, there is only one such free parameter, which controls the degree of locality 
of baryon-antibaryon pair formation. 

In this paper we present an analysis of inclusive charged particle and identified 7r ± , K ± , 
and p/p production in hadronic Z° decays collected by the SLC Large Detector (SLD). The 
analysis is based upon the approximately 400,000 hadronic events obtained in runs of the 
SLC between 1996 and 1998, and supersedes our previous measurements [12]. We measure 
differential production rates for these hadron species in a flavor-inclusive sample of hadronic 
Z° decays and use the results to test the predictions of MLLA QCD+LPHD and of the 
three models just described. We also measure these differential production rates separately 
in Z° decays into light flavors (mm, dd and ss), cc and bb, which provide improved tests of the 
QCD predictions and new tests of the models that separate the heavy hadron production and 
decay modelling from that of the rest of the hadronization process. In addition we measure 
hadron and antihadron differential production rates in light quark jets, thereby obtaining 
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precise information on leading particles and additional stringent tests of the models. 

In section 2 we describe the SLD, in particular its precision vertexing and particle iden- 
tification capabilities. Section 3 describes the selection of hadronic events, their separation 
into samples of different primary flavors using displaced vertices from heavy hadron decays, 
and the selection of light quark and antiquark hemispheres using the large production asym- 
metry in polar angle induced by the polarization of the SLC electron beam. In section 4 we 
present a measurement of the inclusive stable charged particle production rate. We describe 
the hadron species identification and present results for flavor-inclusive events in section 5, 
and section 6 presents results separately for light- (Z° — >uu, dd, ss), c- (Z° ^cc) and 6-flavor 
(Z° — > bb) events. In sections 7 and 8 we use these results to test the predictions of MLLA 
QCD+LPHD, and extract total yields of each hadron species per event of each flavor, re- 
spectively. We present measurements of leading particle production in light-flavor jets in 
section 9, and summarize the results in section 10. 

2 The SLD 

This analysis of data from the SLD [17] used charged tracks measured in the Central Drift 
Chamber (CDC) [18] and upgraded CCD Vertex Detector (VXD3) [19], and identified in the 
Cherenkov Ring Imaging Detector (CRID) [12, 20]. Energy deposits reconstructed in the 
Liquid Argon Calorimeter (LAC) [21] were used in the initial hadronic event selection and in 
the calculation of the event thrust [22] axis. Momentum measurement was provided by an 
axial magnetic field of 0.6T. The CDC and VXD3 gave a momentum resolution of cr p± /p± 
= 0.01 © 0.0026p^, where p± is the track momentum transverse to the beam axis in GeV/c. 

In the plane normal to the beamline the centroid of the micron-sized SLC interaction 
region was reconstructed from tracks in sets of approximately thirty sequential hadronic Z° 
decays to a precision of a r jp ~3 /im and used as an estimate of the primary interaction point 
(IP). The IP position along the beam axis was determined event by event using charged 
tracks, with an average resolution of af P ~ 20 /im. Including the uncertainty on the IP posi- 
tion, the resolution on the charged track impact parameter S was parametrized in the plane 
perpendicular to the beamline as a r ^ =8©29/(psin 3 / 2 9) /im, and in any plane containing 
the beam axis as erf =9©29/G»sin 3/2 #) /im, where 6 is the track polar angle with respect to 
the beamline. 

The barrel CRID covered the polar angle range | cos#| < 0.68, and comprised liquid and 
gaseous radiator systems with refractive indices of 1.278 and 1.00176, respectively. In the 
high momentum limit, an average of 13 (10) Cherenkov photons was reconstructed per track 
traversing the liquid (gas) system, with an average Cherenkov angle resolution of 15 (4.3) 
mrad. The combination of these systems gave efficient and pure separation of pions, kaons 
and protons over much of the kinematic range, as detailed in section 5. 

3 Event Selection 

The trigger and initial selection of hadronic events are described in [23]. The analysis 
presented here is based on charged tracks measured in the CDC and VXD3. A set of cuts 
was applied in order to select charged tracks and events well-contained within the detector 
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uu, dd, ss 
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bb 


uu, dd, ss 


cc 


bb 


lids-tag 


0.703 


0.190 


0.010 


0.925 


0.070 


0.005 


c-tag 


0.061 


0.551 


0.105 


0.204 


0.638 


0.158 


6-tag 


0.001 


0.024 


0.815 


0.005 


0.023 


0.973 


no-tag 


0.235 


0.235 


0.070 


0.721 


0.201 


0.078 



Table 1: Efficiencies for simulated events in the three flavor categories to be tagged as light 
(uds), c or b flavor, or none of these. The three rightmost columns indicate the composition 
of each simulated tagged sample assuming Standard Model relative flavor production. 



acceptance. Tracks were required to have: (i) a closest approach to the beam axis within 
5 cm, and within 10 cm along the beam axis of the measured IP; (ii) a polar angle 9 with 
respect to the beam axis with | cos6>| < 0.80; (iii) a momentum transverse to this axis p± > 
150 MeV/c; and (iv) a momentum p < 50 GeV/c. Events were required to have: a minimum 
of seven such tracks; a visible energy E vis > 18 GeV, calculated from the accepted tracks, 
assigned the charged pion mass; a thrust axis polar angle t with respect to the beam axis 
with | cos# t | < 0.71; and the CDC and VXD3 active and a well-measured IP position. A 
sample of 284,494 events passed these cuts. For the identified hadrons, the CRID was also 
required to be active, giving a sample of 232,802 events. The non-hadronic background was 
estimated to be 0.1%, dominated by Z° — > t + t~ events. 

Samples of events enriched in light, c, and b primary flavors were selected [24] using 
tracks with well measured impact parameters. For each event we defined n sig as the number 
of tracks with an impact parameter greater than three times its estimated error, 5 > 3(7$. We 
also ran a topological vertex finding algorithm [25] on the set of tracks in each hemisphere and 
considered the pt-corrected mass [25] M pt of any secondary vertex (significantly separated 
from the IP) found. Any event containing a vertex with M pt > 2 GeV/c 2 was assigned to 
the 6-tagged sample. If not 6-tagged, an event was c-tagged if either: it contained a vertex 
with 0.5<M pt <2 GeV/c 2 , a total momentum of tracks assigned to the vertex P vtx >2 GeV/c 
and cP v tx — 14:C 2 M pt >— 10 GeV; or it contained no secondary vertex and had 2 < n S i g < 3. 
Events with no secondary vertex and n S i g = were assigned to the light-tagged sample. 
The remaining events were kept as an untagged sample. The 15% of the data taken in 1996 
were excluded due to uncertainties in the simulation, and the light-, c- and 6-tagged samples 
comprised roughly 83,000, 28,000 and 33,000 events, respectively. Selection efficiencies were 
calculated from a detailed detector simulation based on a version of JETSET tuned to the 
world's data on inclusive particle production and D and B hadron production and decay [26]; 
efficiencies and sample purities for the above selection are listed in table 1. 

Samples of hemispheres enriched in light-quark and light-antiquark jets were selected 
by exploiting the large electroweak forward-backward production asymmetry with respect 
to the beam direction induced by the high polarization of the SLC electron beam. Here a 
looser light-flavor event tag of n S i g = was applied, giving a simulated uds efficiency of 84% 
and purity of 89%, with the background dominated by c-flavor events. The event thrust 
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axis was used to approximate the initial qq axis and was signed such that its z-component 
was along the electron beam direction, t z > 0. Events in the central region of the detector, 
where the production asymmetry is small, were removed by the requirement \t z \ > 0.15, 
leaving 109,000 events. The quark-tagged hemisphere in events with left- (right-) handed 
electron beam polarization was defined to comprise the set of tracks with positive (negative) 
momentum projection along the signed thrust axis. The remaining tracks in each event 
were defined to be in the antiquark-tagged hemisphere. For the selected event sample, the 
average magnitude of the e~ beam polarization was 0.73. Using this value and assuming 
Standard Model couplings, a tree- level calculation gives a quark (antiquark) purity of 0.73 
in the quark- (antiquark-)tagged sample. 

4 Stable Charged Particle Rates 

We first measured the differential production rate of all stable charged particles defined con- 
ventionally as the sum of electrons, muons, pions, kaons and protons that are either primary 
hadrons or products of a chain of decays of particles with proper lifetime less than 3xl0~ 10 s. 
Tracks satisfying the requirements in sec. 3 were counted in momentum bins and corrected, 
using our detector simulation, for the track finding and selection efficiency and resolution, 
non-hadronic event background, and spurious tracks from beam-related backgrounds and 
interactions in the detector material. The resulting integrated rate was constrained to be 
20.95±0.21 charged particles per event, an average of measurements [10, 27, 28] in Z° decays. 

The momentum dependence of the selection was constrained by comparing the properties 
of reconstructed charged tracks in measured and simulated r lepton decays, for which the 
momentum and particle type distributions are well known. We selected e + e _ — > t + t~ 
events according to the criteria in [29] , with the additional requirement of at least one track 
in the event with p >7 GeV/c, to reduce beam-related and two-photon event background. 
Comparisons were made for 1-, 3-, and 5-prong decays separately, as well as for hemispheres 
in which 2 or 4 tracks were found, giving constraints on both isolated tracks and those in 
close proximity to others. We also checked the momentum distributions for tagged electrons, 
muons, pions, and kaons, and the overall multiplicity distribution in selected r-pair events. 
In all cases the simulation was consistent with the data. Linear fits were made to the ratios 
of data:simulated tracks, and a momentum-dependent uncertainty of O.lOx \p — 2.7 GeV/c|% 
was assigned, reflecting the error on a typical fitted slope. Here p is the particle momentum 
in GeV/c and 2.7 GeV/c is the average momentum of all charged particles, so that this 
uncertainty is independent of the normalization uncertainty noted above. Variations of the 
background levels, detector material and momentum resolution were also considered and 
found to have much smaller effects, except at the lowest and highest momenta. 

The inclusive charged particle differential production rate is listed in table 2, in terms 
of the variables momentum p, scaled momentum x p = 1pjE CM and £ = — \wx p . The x p 
distribution is shown in fig. 1, and compared with the predictions of the JETSET, UCLA 
and HERWIG event generators described in section 1, using the default parameter values 
for each model. The JETSET model is the most consistent with the data, but predicts a 
slightly softer spectrum; the UCLA model describes the data over most of the x p region, 
but falls low for x p >0.7 and rr p <0.015 units; the HERWIG model predicts a spectrum quite 
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different from that seen in the data. Our results for the shape of the spectrum are consistent 
with those published previously [10, 27, 28]. 

5 Charged Pion, Kaon and Proton Production 

Additional track selection cuts [24] were applied to remove tracks that interacted or scattered 
through large angles before exiting the CRID and to ensure that the CRID performance was 
well-modelled by the simulation. Good information from the liquid (gas) system was required 
for tracks in the "liquid" ("gas") region, with momenta below 2.25 GeV/c (above 10.5 
GeV/c); good information from both systems was required for momenta in the "combined" 
region, 2.25-10.5 GeV/c. Tracks were required to have at least 40 CDC hits, at least one 
of which was at a radius of at least 92 cm, to extrapolate through an active region of 
the appropriate radiator (s), and to have at least 40 (70)% of their expected liquid (gas) ring 
contained within a sensitive region of the CRID TPCs. The latter requirement included dead 
and inefficient regions within the TPCs and also rejected tracks with p > 2.25 GeV/c for 
which there was a saturated CRID hit within a 2 cm radius (~the maximum ring radius) of a 
point 1 cm behind the expected gas ring center in the photoelectron drift direction. Saturated 
hits arise from the passage of tracks through the TPCs and from various backgrounds, 
and mask single photoelectron hits in their vicinity. Liquid region tracks that extrapolated 
through an active TPC were required to have a saturated hit within 1 cm of the extrapolated 
track, indicating that the track traversed the TPC and hence the liquid radiator; gas region 
tracks were required to have either such a saturated hit or the presence of at least four hits 
consistent with a liquid ring. These cuts accepted 71%, 68% and 74% of the tracks within 
the CRID geometrical acceptance in the liquid, combined and gas regions, respectively. 

Tracks were identified using a global likelihood technique [12, 30]. For each track and each 
of the five hypotheses i — e, /i, n, K, p, a likelihood Lj was calculated based upon: the number 
of detected photoelectrons and their measured Cherenkov angles; the expected number of 
photons; the expected Cherenkov angle; and terms accounting for random backgrounds and 
hits consistent with Cherenkov radiation from other tracks in the event. The best hypothesis 
for each track was used to determine its contribution to the background for other tracks, 
and the calculation was iterated until there was no change in any best hypothesis. Particle 
separation was based upon differences between logarithms of the three likelihoods, Ci = In Lj, 
% — 7r, K, p. A track in the liquid (combined, gas) region was identified as species j if Cj 
exceeded both of the other log-likelihoods by at least 3 (2) units. The electron and muon 
likelihoods are generally quite similar to the pion likelihood, and the leptons were included 
in the pion category at this stage. 

We quantified the identification performance in terms of a momentum-dependent iden- 
tification efficiency matrix E, each element Eij of which represents the probability that the 
selected track from a true i-hadron is identified as a j-hadron, with i,j = n,K,p. The ele- 
ments of this matrix were determined where possible from the data [24, 31]. For example, 
tracks from selected K Q S and r decays were used as "pion" test samples, having estimated 
kaon plus proton contents of 0.3% and 1.7% respectively. Probabilities for such tracks to be 
identified as pions, kaons and protons were compared with the predictions of our detailed 
MC simulation, which was found to describe the momentum dependence of the efficiencies 
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q c;n 
o.OU— 


n kc\ 

y.ou 


8.982 


0.2687±0 


0011±0 


0032±0 


nnnQ 
UUUo 


12.25±0 


15 


2.328±0.029 


n kc\ 

y.ou— 


1 n c;n 
1U.0U 


9.980 


0.2120±0 


0010±0 


0026±0 


nnno 
UUUz 


9.67±0 


12 


2.039±0.027 


in c 

1U.0— 


11.0 


10.99 


0.1700±0 


0009±0 


0022±0 


nnno 
UUUz 


7.75±0 


11 


1.792±0.025 


11.0— 


It) c 

1Z.0 


11.99 


0.1351±0 


000^ 


3±0 


001^ 


3±0 


nnm 
UUU1 


6.161±0 


088 


1.551±0.023 


1 o c: 
lz.0— 


lo.O 


12.99 


0.1103±0 


0007±0 


0016±0 


nnm 
UUU1 


5.029±0 


076 


1.368±0.022 


lo.O— 


14.0 


13.99 


0.0889±0 


0006±0 


0013±0 


nnm 
UUU1 


4.053±0 


065 


1.184±0.020 


1 A K 
14.0— 


1 n 
10. U 


15.22 


0.0688±0 


0005±0 


0011±0 


nnm 
UUU1 


3.139±0 


052 


0.997±0.017 


16.0- 


17.5 


16.72 


0.0513±0 


0004±0 


0009±0 


0001 


2.338±0 


042 


0.812±0.016 


17.5- 


19.0 


18.23 


0.0383±0 


0003±0 


0007±0 


0001 


1.748±0 


034 


0.661±0.014 


19.0- 


20.5 


19.72 


0.0291±0 


0003±0 


0006±0 


0001 


1.326±0 


028 


0.544±0.012 


20.5- 


22.0 


21.21 


0.0221±0 


0003±0 


0005±0 


0001 


1.008±0 


023 


0.444±0.011 


22.0- 


24.0 


22.96 


0.0159±0 


0002±0 


0004±0 


0001 


0.724±0 


018 


0.3468±.0090 


24.0- 


26.0 


24.94 


0.0105±0 


0001±0 


0003±0 


0001 


0.480±0 


013 


0.2513±.0073 


26.0- 


30.0 


27.74 


0.0063±0 


0001±0 


0002±0 


0001 


0.285±0 


009 


0.1700±.0056 


30.0- 


35.0 


32.14 


0.0025±0 


0000±0 


0001±0 


0001 


0.114±0 


005 


0.0800±.0035 


35.0- 


45.0 


38.23 


0.0005±0 


0000±0 


0000±0 


0001 


0.024±0 


001 


0.0207±.0013 



Table 2: Production rate for all stable charged particles in terms of momentum p, scaled 
momentum x p and £ = — lnx p . For momentum we give statistical errors and the systematics 
arising from track-finding efficiency (including an overall 1% normalization term) and the 
sum of other sources, which is dominated by backgrounds (momentum resolution) at low 
(high) momentum; in the other columns these have been summed in quadrature. 
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■ Systematic Error Band f 
UCLA:JETSET / 




Figure 1: a) Charged track x p distribution in hadronic events (dots) with logarithmic and 
(inset) linear vertical scales, compared with the predictions of three models (lines), b) 
The same data divided by the JETSET prediction. The error bars in (b) are statistical 
and the shaded band represents the systematic uncertainty; all errors except an overall 1% 
normalization uncertainty are included in (a). 
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well and to reproduce their amplitudes to within a few percent. Functional forms, chosen to 
describe the momentum dependence of both data and simulated test samples as well as that 
of simulated true pions in hadronic events, were fitted to the data, except for momenta below 
0.8 GeV/c where there was significant structure on the scale of the bin size and bin-by-bin 
corrections were used. The simulation was used to correct the fitted parameters for non-pion 
content in the test samples and differences in tracking performance between tracks in these 
samples and those passing our selection cuts in hadronic events. The resulting identification 
efficiency functions, E nn , E n K and E np , are shown in the leftmost column of fig. 2. A similar 
procedure using only 7r and p likelihoods was used to measure the 7r-p separation in the 
liquid (gas) system for p > 2 (17) GeV/c. This information was combined with that from 
test samples of protons and kaons from decays of A hyperons and mesons, respectively, 
to derive the remaining efficiencies in fig. 2. 

The bands in fig. 2 encompass the upper and lower systematic error bounds on the 
efficiencies. There are discontinuities between the liquid, combined and gas regions, and 
there are strong point-to-point correlations within each region. For the diagonal elements, 
these errors correspond to those on the fitted parameters, and fall into four categories. The 
overall amplitude of the efficiency is driven by the average and rms of the number of photons 
detected per track; the corresponding uncertainty is common to all momenta in a given region 
and is also correlated between the liquid and combined regions for protons and the combined 
and gas regions for pions. The Cherenkov angle resolution affects the positions of the falling 
edges; its uncertainty is therefore correlated across the range of a given edge, as well as 
between pions and kaons in the regions 1.5-2.25 GeV/c and 12-45 GeV/c, and between kaons 
and protons in the region 3-8 GeV/c. Performance near a Cherenkov threshold depends on 
the relevant index of refraction and its stability; since the efficiencies change rapidly on the 
scale of our bin sizes, bin-by-bin calibrations were done for E^k (E pp ) in the range p <1 
GeV/c (p < 1.5 GeV/c), and Ekk and E pp in the range 7.5< p <10.5 GeV/c, which are 
completely independent. In the region 10.5-18 GeV/c, protons are below threshold in the 
gas, whereas pions and kaons are well above threshold; here E pp depends largely on the 
background level, and its uncertainty is correlated across this region. 

For the off-diagonal elements, representing misidentification rates, the errors on fitted 
parameters were also used, but subject to a minimum value of 0.0025 to account for the 
limited statistics of the test sample constraints on the momentum dependences. These un- 
certainties correspond to a combination of the effects listed above for the diagonal elements, 
and each is typically dominated by one effect in a given momentum region. These errors 
should also therefore be considered strongly positively correlated across each of the liquid, 
combined and gas regions. 

The identification efficiencies in fig. 2 peak near or above 0.9 and are greater than 0.8 
over wide ranges. The pion (kaon and proton) coverage is continuous from 0.25 GeV/c (0.65 
GeV/c) up to the beam energy, although the efficiencies fall below 0.2 for pions and kaons 
above about 30 GeV/c, and kaons and protons in the 6-9 GeV/c range. Misidentification 
rates are typically at the few percent level, with peak values of up to 0.1. 

In each momentum bin we measured the fractions of the selected tracks that were identi- 
fied as pions, kaons and protons. The observed fractions were related to the true production 
fractions by an efficiency matrix, composed of the values shown in fig. 2. This matrix was 
inverted and used to unfold our observed identified hadron fractions. This analysis procedure 
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Figure 2: Calibrated identification efficiencies for selected tracks. The half- widths of the 
grey bands represent the systematic uncertainties, which are strongly correlated between 
momenta. The off-diagonal efficiencies have been scaled by factors of five for clarity. The 
white bands in the four upper left plots represent the 2x2 matrix used below kaon threshold 
in the liquid system. 
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does not require that the sum of the charged hadron fractions be unity; instead the sum was 
used as a consistency check, which was found to be satisfied at all momenta (see fig. 3). For 
momenta below 0.65 GeV/c, we could not distinguish kaons from protons, but pions could 
be identified down to 0.25 GeV/c. An analogous 2x2 analysis of pions and non-pions was 
used in this region, and we present only the pion fraction. 

The background from electrons and muons was estimated from the simulation to be about 
5% of the tracks in the inclusive flavor sample, predominantly from c- and 6-flavor events. 
The fractions were corrected for the lepton backgrounds using the simulation, as well as for 
the effects of beam-related backgrounds, particles interacting in the detector material, and 
particles decaying outside the tracking volume. The conventional definition of a final-state 
charged hadron was used, namely a charged pion, kaon or proton that is either from the 
primary interaction or a product of a chain of decays of hadrons each with a proper lifetime 
less than 3xl0~ 10 s. 

The measured charged hadron fractions in inclusive hadronic Z° decays are shown in fig. 3 
and listed in tables 3-5. The systematic errors were determined by propagating the errors 
on the calibrated efficiency matrix and are strongly correlated as described above. They 
are indicated by the shaded regions in fig. 3. The errors on the points below ~15 GeV/c 
are dominated by the systematic uncertainties; those above ~30 GeV/c are dominated by 
statistical uncertainties. The sum of fractions is consistent with unity everywhere. 

Pions are seen to dominate the charged hadron production at low momentum, and to 
decline steadily in fraction as momentum increases. The kaon fraction rises steadily, ap- 
proaching the pion fraction at high momentum. The proton fraction rises to a plateau value 
of about one-tenth at about 10 GeV/c, then declines at the highest momenta. 

In fig. 4 we compare our results with the predictions of the JETSET, UCLA and HERWIG 
models. The momentum dependence for each of the three hadron species is reproduced 
qualitatively by all models. All three models overestimate the kaon fraction for momenta 
below about 1.5 GeV/c, and underestimate it for momenta above about 15 GeV/c. The 
UCLA model overestimates the pion fraction by about 2a (taking into account the correlation 
in the experimental errors) in the 5-15 GeV/c range, but is the only model consistent with 
the behavior above about 20 GeV/c. The JETSET model overestimates the proton fraction 
at all momenta, but describes the momentum dependence; HERWIG and UCLA predict 
behavior at high momentum that is inconsistent with the proton data. 

Since the JETSET model has a number of parameters that control specific aspects of 
relative particle production, we have investigated some simple changes. We find that reduc- 
ing the diquark:quark ratio (PARJ(l)) from 0.10 to 0.08, while leaving all other parameters 
unchanged, results in a good description of the proton fraction. The kaon fraction is sen- 
sitive to both the s-quark probability (PAR J (2)) and the strange vector :pseudoscalar ratio 
(PARJ(12)), however no combination of these two parameters was found to give a good 
description of the kaon fraction over the entire measured momentum range. 

In order to obtain charged hadron production rates, the simulation was used to subtract 
the contribution of all particles (mostly leptons) that were not charged pions, kaons or 
protons from our measured total charged production rate (table 2). The pion, kaon and 
proton fractions were multiplied by this adjusted rate to obtain the individual rates tabulated 
in tables 3-5 as a function of scaled momentum. In fig. 5 we compare our £ distributions with 
the predictions of the three models. The features for pions are similar to those for all charged 
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Figure 3: a) Measured charged hadron production fractions in hadronic Z° decays. The 
circles represent the 7r =l= fraction, the squares the K ± fraction, the diamonds the p/p fraction. 
The error bars are statistical only; the shaded areas indicate the systematic errors, and are 
connected across momentum regions where there is a strong positive correlation, b) The 
sum of the three fractions. 
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x p Range 


<x p > 


U 


(1/N evts ) dn^/dxp 





005- 


-0 


008 





0066 





951±0 


002±0 


010 


471. 8±1 


3 


±9 


2 





008- 


-0 


010 





0088 





933±0 


001±0 


007 


470.4±1 


1 


±6 


6 





010- 





012 





0109 





910±0 


001±0 


006 


434.6±1 


1 


±5 








012- 


-0 


014 





0131 





901±0 


001±0 


006 


388.8±1 





±4 








014- 


-0 


016 





0153 





885±0 


001±0 


006 


352.7±0 


9 


±3 


3 





016- 


-0 


022 





0191 





888±0 


001±0 


004 


294.8±0 


5 


±2 


2 





022- 


-0 


027 





0246 





870±0 


001±0 


003 


229.6±0 


5 


±1 


3 





027- 


-0 


033 





0301 





860±0 


001±0 


003 


185.0±0 


4 


±0 


9 





033- 





038 





0356 





845±0 


001±0 


003 


150.6±0 


4 


±0 


7 





038- 


-0 


044 





0411 





831±0 


002±0 


005 


125. 6±0 


4 


±0 


9 





044- 


-0 


049 





0465 





828±0 


002±0 


008 


106.5±0 


4 


±1 


1 





049- 


-0 


055 





0521 





810±0 


002±0 


012 


90.40±0 


35 


±1 


33 





055- 


-0 


060 





0576 





794±0 


002±0 


009 


77.38±0 


31 


±0 


91 





060- 


-0 


066 





0630 





784±0 


002±0 


008 


67.39±0 


29 


±0 


70 





066- 


-0 


071 





0685 





780±0 


002±0 


010 


59.40±0 


27 


±0 


75 





071- 


-0 


077 





0740 





768±0 


003±0 


009 


52.57±0 


25 


±0 


60 





077- 


-0 


082 





0795 





760±0 


003±0 


008 


46.76±0 


24 


±0 


50 





082- 


-0 


088 





0850 





752±0 


003±0 


008 


41.70±0 


23 


±0 


43 





088- 


-0 


099 





0931 





738±0 


002±0 


007 


35.26±0 


15 


±0 


36 





099- 


-0 


110 





1040 





726±0 


002±0 


007 


28.89±0 


13 


±0 


29 





110- 


-0 


121 





1150 





712±0 


003±0 


007 


23.88±0 


12 


±0 


25 





121- 


-0 


143 





1310 





702±0 


002±0 


007 


18.69±0 


08 


±0 


19 





143- 


-0 


164 





1530 





696±0 


002±0 


006 


13.85±0 


07 


±0 


14 





164- 


-0 


186 





1750 





673±0 


003±0 


006 


10.16±0 


06 


±0 


11 





186- 


-0 


208 





1970 





662±0 


003±0 


004 


7.812±0 


050 


±0 


069 





208- 


-0 


230 





2189 





653±0 


004±0 


004 


6.076±0 


044 


±0 


061 





230- 


-0 


252 





2410 





629±0 


004±0 


005 


4.674±0 


039 


±0 


053 





252- 


-0 


274 





2629 





616±0 


005±0 


005 


3.632±0 


035 


±0 


044 





274- 


-0 


296 





2849 





601±0 


005±0 


004 


2.886±0 


031 


±0 


037 





296- 


-0 


318 





3068 





594±0 


006±0 


004 


2.292±0 


028 


±0 


031 





318- 


-0 


351 





3338 





586±0 


006±0 


009 


1.749±0 


021 


±0 


034 





351- 


-0 


384 





3666 





577±0 


007±0 


010 


1.275±0 


018 


±0 


028 





384- 


-0 


417 





3997 





557±0 


008±0 


010 


0.921±0 


016 


±0 


022 


u 


A 1 1 

41 (- 


-U 




u 




u 


04z±U 


U1U±U 


nin 

U1U 


U.boU±U 


014 


±0 


Ulo 





450- 


-0 


482 





4651 





523±0 


012±0 


Oil 


0.499±0 


013 


±0 


014 





482- 


-0 


526 





5035 





491±0 


013±0 


010 


0.338±0 


010 


±0 


010 





526- 


-0 


570 





5470 





492±0 


018±0 


Oil 


0.226±0 


009 


±0 


007 





570- 


-0 


658 





6083 





469±0 


018±0 


012 


0.130±0 


005 


±0 


005 





658- 


-0 


768 





7047 





467±0 


032±0 


018 


0.0526±0 


0037±0 


0029 





768- 


-1 


000 





8383 





479±0 


074±0 


048 


0.0113±0 


001^ 


3±0 


0013 


0.005- 


-1.000 




15.74±0.01 


±0.17 



Table 3: Charged pion fraction /V and differential production rate (1/N e ^ s )dn 7r / 'dx p per 
hadronic Z° decay. The first error is statistical, the second systematic. <x p > is the average 
rrp-value of charged tracks in each bin. The last row gives the integral over the x p range of 
the measurement. A 1.0% normalization uncertainty is included in the systematic error on 
the integral, but not in those on the differential rate. 
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Xp Range 


<x p > 


h 


(i/rw 


) dn K /dx p 





1)14- 





016 


A 




n 1 c o 

0153 


A 




nTo i n 

072±0 


nno_i_n 

002±0 


noo 

023 


no rn i r\ 


r a 


I A 

±9 


0£? 

26 


A 




A1 

Olb- 


A 

-0 


022 


A 




0191 


A 




r\n r i a 

065±0 


An 1 i n 

001±0 


nn r 

005 


zl.o ^ ±U 


on 

20 


1 1 

±1 


PT T 

57 


A 




Uzz- 


A 

-0 


noT 

027 


A 




no a p 


A 




noo_i_n 

082±0 


nn i i n 

001±0 


nno 

003 


Zl.0Z±U 


i n 

19 


I A 

±0 


on 

80 


A 




A O T 

027- 


-0 


033 


A 




n on 1 

0301 


A 




A A 1 I A 

091±0 


nn 1 i n 

001±0 


nno 

002 


iy.oo±u 


1 O 

18 


i n 

±0 


r o 

53 


A 




A O O 

033- 


A 

-0 


A o o 

038 


A 




n o p* r* 
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A 




1 A 1 1 A 

101±0 


A A 1 I A 

00l±0 


nno 

002 


1 O no i n 

18.02±0 


1 Z" 1 

16 


1 A 

±0 


A A 

44 


A 




A O O 

038- 


A 

-0 


1 \ A A 

044 


A 




n a i i 

0411 


A 




1 1 A I A 

114±0 


A A 1 1 A 

001±0 


nno 

003 


1 1 .1 1 ±0 


1 T 

17 


1 A 

±0 


A O 

43 


a 




1 \ A A 

044- 


A 

-0 


A A A 

049 


A 




a 1 n r 

0465 


A 




1 OO I A 

123±0 


A A 1 I A 

001±0 


nn a 

004 


15. /8±0 


1 T 

17 


1 A 

±0 


A T 

47 


A 


A 1 A 

049- 


A 

-0 


055 


A 




0521 


A 




1 01 I A 

131±0 


nno i n 

002±0 


nn a 

004 


1 /! c^/i i n 

14.oo4±0 


i n /i 
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i n 

±0 


A A O 

442 


a 

u 


nrr 

055- 


A 

-0 


A A 

060 


A 




n r Tr 

0576 


A 




i on i n 

139±0 


nno i n 

002±0 


nn r 

005 




i on 

189 


i n 

±0 


r n o 

503 


A 




A A 

060- 


A 

-0 


A/" 1 /" 1 

066 


A 




n r* o n 

0630 


A 




1 A T 1 A 

147±0 


nno i n 

002±0 


A A/"" 1 

006 


1 O rrnn i n 

12.o99±0 


1 T/" 1 
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1 A 

±0 


V V o 

558 


a 

o 


066- 


A 

-0 


071 


A 




n/"* o r 

0685 


A 




1 r o i n 

158±0 


nno i n 

002±0 


nno 

008 


1 O nif i n 

12.03o±0 


165 


i n 

±0 


635 


A 

(J 


071- 


A 

-0 


077 


A 




nT /i n 

0740 


A 




1 CP i n 

166±0 


nno i n 

002±0 


nnn 

009 


1 1 Q /in i n 


1 /"* O 
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i n 

±0 


/^OO 
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A 
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077- 


-0 


A O O 
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A 




n Tn r 

0795 


A 




1 f f 1 a 

166±0 


nno i n 

003±0 


n i n 

010 


1U.ZU 1 ±U 


1 P A 
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i n 

±0 


/ ■ A ' ) 
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A 


A O O 

082- 


-0 


A O O 

088 


A 




A o r A 

0850 


A 




i to i n 

172±0 


nno i n 

003±0 


n i n 

010 


9.0 ( 1±U 


160 


i n 

±0 
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A 




A O O 

088- 


-0 


AAA 

099 


A 




nno 1 

0931 


A 




1 o 1 in 

181±0 


nno i n 

002±0 


n 1 1 

011 


o.o / 1±U 


1 1 o 

113 


i n 

±0 


r n r 

505 


A 




AAA 

099- 


A 

-0 


1 1 A 
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A 




1 n a n 

1040 


A 




1 A {' 1 A 

196±0 


nno i n 

003±0 


n 1 1 

011 




11/1 

114 


1 A 

±0 


/i A n 

440 


A 




1 1 A 

110- 


A 

-0 


1 O 1 

121 


A 




1150 


A 




1 1 A 

216±0 


nn /i i n 

004±0 


n I o 

012 


l .2,6 1 ±0 


i on 

120 


1 A 

±0 


OA P" 

395 


A 




121- 


A 




143 


A 




1 o i n 

1310 


A 




216±0 


nno i n 

003±0 


n 1 a 

014 


5. /4o±0 


non 

089 


i n 

±0 


o/^n 

369 


a 

u 


143- 


A 




1 {> A 

164 


A 




1 r o n 

1530 


A 




1 A A 1 A 

199±0 


nn r i n 

005±0 


n i n 

019 


Q nrrn i n 


i no 

102 


i n 

±0 


O O 1 

381 


A 




1 {' A 

164- 


A 

-0 


i o r 1 

186 


A 




1750 


A 




OO A 1 A 

230±0 


nnn i n 

009±0 


A O V 

035 
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o.4/o±U 


1 O A 

134 


1 A 

±0 


V o o 

532 


A 




1 O/" 1 

186- 


A 

-0 


O AO 

208 


A 




1 ATA 

1970 


A 




OOO 1 A 

232±0 


A A /I 1 A 

004±0 


A O V 

035 


o 7onj_n 
2. M9±0 


n /i t 

047 


1 A 

±0 


/i i n 

419 
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208- 


A 




230 


A 




o i o n 

2189 


A 




or* a i n 

264±0 


nn a i n 

004±0 


AIT 

017 


O A rro i n 

z.4oz±U 


no t 

037 


i n 

±0 


1 /^ O 

163 


A 




230- 


A 

-0 


252 


A 




o /i i n 

2410 


A 




o n c in 

256±0 


nn /i i n 

004±0 


nno 

008 


1 nno i n 
1.903±0 


non 

030 


i n 

±0 


C\P O 

063 





252- 


A 

-0 


274 


A 




o/^on 

2629 


A 




A/ - t i n 

267±0 


nn /i i n 

004±0 


006 


1 c7i in 

1.0 / 4±U 


no t 

027 


i n 

±0 


no /"* 

036 


A 




O T /I 

274- 


A 




296 


A 




O O .1 A 

2849 


A 




ooo i n 

283±0 


nn r i n 

005±0 


nn /i 

004 


1 Q/^n i n 

l.ooU±U 


no a 

024 


i n 

±0 


nor 

026 


A 




a A z" 1 

296- 


A 

-0 


o i o 

318 


A 




O A/"" 1 O 

3068 


A 




OA A 1 A 

290±0 


A A P* - 1 A 

005±0 


nn /i 

004 


1 11 o i n 

l.llo±U 


noo 

022 


1 A 

±0 


AO A 

020 


A 




318- 


A 

-0 


351 


A 




o o o o 

3338 


A 




OA O 1 A 

298±0 


aa r 1 A 

005±0 


nn a 

004 


n onnin 
0.890±0 


016 


1 A 

±0 


A 1 T 

017 




on 

351- 


-0 


384 





o/" 1 /" 1 c 

3666 





onn i n 

309±0 


AAP 1 A 

006±0 


006 


O.o83±0 


014 


±0 


A1 /" 

016 





384- 


-0 


417 





3997 





343±0 


007±0 


008 


U.Ob ( ±U 


013 


±0 


015 





417- 


-0 


450 





4325 





345±0 


008±0 


009 


0.433±0 


012 


±0 


014 





450- 


-0 


482 





4651 





368±0 


010±0 


011 


0.351±0 


011 


±0 


012 





482- 


-0 


526 





5035 





384±0 


011±0 


013 


0.264±0 


008 


±0 


010 





526- 


-0 


570 





5470 





411±0 


016±0 


015 


0.188±0 


008 


±0 


008 





570- 


-0 


658 





6083 





439±0 


017±0 


017 


0.122±0 


005 


±0 


006 





658- 


-0 


768 





7047 





431±0 


032±0 


018 


0.0485±0 


0037±0 


0027 





768- 


-1 


000 





8383 





328±0 


090±0 


042 


0.0078±0 


0022±0 


0011 


0.014- 


-1.000 




2.074±0.006 


±0.066 



Table 4: Charged kaon fraction and differential production rate per hadronic Z° decay. 
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Xp Range 


<x p > 


/p 


(1/N e „ ts ) dn p /dx p 





014- 


-0 


016 





0153 





036±0 


001±0 


013 


14.51±0 


52 


±5 


08 





016- 


-0 


022 





0191 





052±0 


001±0 


008 


17.32±0 


27 


±2 


58 





022- 





027 





0246 





052±0 


001±0 


009 


13.75±0 


29 


±2 


50 





027- 


-0 


033 





0301 





052±0 


001±0 


006 


11.12±0 


17 


±1 


24 





033- 


-0 


038 





0356 





060±0 


001±0 


003 


10.75±0 


14 


±0 


47 





038- 


-0 


044 





0411 





060±0 


001±0 


002 


9.048±0 


123 


±0 


350 





044- 


-0 


049 





0465 





060±0 


001±0 


002 


7.669±0 


117 


±0 


298 





049- 


-0 


055 





0521 





066±0 


001±0 


003 


7.410±0 


113 


±0 


294 





055- 


-0 


060 





0576 





068±0 


001±0 


003 


6.587±0 


109 


±0 


259 





060- 


-0 


066 





0630 





067±0 


001±0 


003 


5.788±0 


105 


±0 


238 





066- 


-0 


071 





0685 





070±0 


001±0 


003 


5.344±0 


100 


±0 


228 





071- 


-0 


077 





0740 





073±0 


002±0 


003 


4.987±0 


104 


±0 


229 





077- 


-0 


082 





0795 





069±0 


002±0 


004 


4.278±0 


100 


±0 


242 





082- 


-0 


088 





0850 





074±0 


002±0 


005 


4.117±0 


101 


±0 


253 





088- 


-0 


099 





0931 





076±0 


002±0 


006 


3.633±0 


072 


±0 


269 





099- 


-0 


110 





1040 





076±0 


002±0 


008 


3.036±0 


076 


±0 


300 





110 


-0 


121 





1150 





077±0 


002±0 


Oil 


2.568±0 


081 


±0 


357 





121- 


-0 


143 





1310 





081±0 


003±0 


015 


2.165±0 


069 


±0 


398 





143- 


-0 


164 





1530 





097±0 


005±0 


023 


1.931±0 


096 


±0 


452 





164- 





186 





1750 





106±0 


009±0 


039 


1.603±0 


133 


±0 


594 





186- 


-0 


208 





1970 





074±0 


004±0 


022 


0.871±0 


045 


±0 


255 





208- 





230 





2189 





098±0 


003±0 


019 


0.912±0 


030 


±0 


179 





230- 


-0 


252 





2410 





104±0 


003±0 


008 


0.775±0 


025 


±0 


062 





252- 





274 





2629 





108±0 


004±0 


007 


0.639±0 


022 


±0 


044 





274- 


-0 


296 





2849 





106±0 


004±0 


007 


0.511±0 


019 


±0 


033 





296- 


-0 


318 





3068 





109±0 


004±0 


006 


0.419±0 


016 


±0 


024 





318- 


-0 


351 





3338 





120±0 


004±0 


006 


0.358±0 


Oil 


±0 


018 





351- 


-0 


384 





3666 





115±0 


004±0 


005 


0.254±0 


009 


±0 


012 





384- 


-0 


417 





3997 





105±0 


004±0 


005 


0.173±0 


008 


±0 


008 





417- 


-0 


450 





4325 





113±0 


005±0 


004 


0.141±0 


007 


±0 


005 





450- 


-0 


482 





4651 





100±0 


006±0 


003 


0.0950±0 


0055±0 


0036 





482- 


-0 


526 





5035 





100±0 


005±0 


003 


0.0688±0 


0039±0 


0027 





526- 


-0 


570 





5470 





103±0 


007±0 


003 


0.0470±0 


0032±0 


0018 





570- 


-0 


658 





6083 





087±0 


006±0 


003 


0.0241±0 


0017±0 


0012 





658- 


-0 


768 





7047 





083±0 


009±0 


004 


0.0093±0 


0010±0 


0006 





768- 


-1 


000 





8383 





062±0 


014±0 


005 


0.0015±0 


0003±0 


0001 


0.014- 


-1.000 




0.984±0.006 


±0.035 



Table 5: Proton plus antiproton fraction and differential production rate per hadronic Z' 
decay. 
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8675A3 Momentum (GeV/c) 

Figure 4: Comparison of our measured charged hadron fractions (symbols) with the pre- 
dictions of the JETSET (dashed lines), UCLA (solid lines) and HERWIG (dotted lines) 
models. 
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1 2 3 4 5 



S ^=-ln(xp) 

Figure 5: Comparison of our measured (a) charged particle and pion and (b) kaon and proton 
production rates (symbols) with the predictions of the JETSET (dashed lines), UCLA (solid 
lines) and HERWIG (dotted lines) models. 
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tracks: all models describe the data qualitatively; JETSET is within about 2% except for 
£ < 1.5 and £ > 4.5; UCLA shows a spectrum very similar to the data but shifted slightly 
toward lower £ values; HERWIG is several percent high (low) for £ < 1 and 3.3 < £ < 4.3 
(1.5 < £ < 2.5). For the protons we see features similar to those seen in the proton fraction 
(fig. 4): JETSET describes the shape of the spectrum but is high by about 20%; UCLA 
describes the data well except for the structure for £ < 0.8; HERWIG is consistent with the 
data for £ > 1.5, but shows a pronounced structure at lower values of £ that is inconsistent 
with the data. For kaons, all three predictions are too high for £ > 3.3, especially HERWIG, 
which is also quite high for £ < 0.6 and low for 1.2 < £ < 2.4; the other two predictions 
are low for £ < 0.8, and consistent with the data elsewhere, although they peak at higher £ 
values than the data. 

Our fractions and production rates are generally consistent with those from previous 
experiments at the Z° [27, 28, 32]. The ALEPH pion spectrum lies above ours at large £; 
the OPAL proton spectrum is lower than all others at small £; all other differences are within 
two standard deviations, considering correlations in the systematic errors. Our measurement 
is the most precise in several regions, most notably for pions in the range 3 < £ < 4, kaons in 
the range 2.5 < £ < 4 and protons in the range 2.5 < £ < 3.5. Measurements based on ring 
imaging and those based on ionization energy loss rates cover complementary momentum 
ranges and could be combined to provide continuous, high-precision coverage over the range 
from the beam momentum (£ = 0) down to p — 0.22 GeV/c (£ = 5.3). 



6 Flavor-Dependent Analysis 

The analyses described above were repeated on the light-, c-, b-, and untagged event samples 
described in section 3, to yield differential production rates i?jf a9 , k = I, c, b, un, respectively; 
the correction for leptons was not performed at this point, so that h = e/iir, K, p. True 
differential rates R™ in events of the three flavor types, m — I, c, b, representing events of 
the types Z° — > uu, dd, ss, Z° — > cc, and Z° — > bb, respectively, were extracted by solving 
for each species h the relations: 

r>ktag mk € mk F m R h 

Rh = ^ m e mk F m ■ (1) 

Here, F m is the fraction of hadronic Z° decays of flavor type m, taken from the Standard 
Model, e mk is an element of the event tagging efficiency matrix (see table 1), and 
represents the momentum-dependent bias of tag k toward selecting events of flavor m that 
contain hadrons of species h. Ideally all biases would be unity in this formulation. The 
biases were calculated from the simulation as B h mk = (n^ tktag / N mtktag ) / (n^/ N m ) , where N m 
(n^J is the number of simulated events (hadrons of species h in events) of true flavor m and 
N m ,ktag ( n m,ktag) ^ s ^ ne num ber of (/i-hadrons in) those events that are tagged as flavor k. The 
diagonal bias values [24] are within a few percent of unity reflecting the small dependences of 
the flavor tags on the charged multiplicity of the event and tracks from strange mesons and 
baryons that decay close to the IP. The sum of the products of biases and efficiencies must 
be unity for events of a given true flavor m, Y> k Bmk e mk = 1; some off-diagonal bias values 
therefore deviate substantially from unity since the corresponding mistagging rate can be 
small, but these do not affect the results as the figure of merit is the diagonal element. 



19 



The unfolded pion rates were corrected for the contributions from leptons by subtracting 
the absolute lepton rate predicted by the simulation. At low momentum this background 
is dominated by electrons from photon conversions, and is a 6% contribution at 1 GeV/c 
that falls rapidly with increasing momentum. For the heavy flavors, electrons and muons 
from semi-leptonic decays of heavy hadrons cause the correction to increase with momentum 
above about 5 GeV/c, reaching 13% (60%) for c-(6-)flavor events in the highest momentum 
bins. 

The resulting differential production rates are listed in tables 6-8. The systematic errors 
listed are only those relevant for the comparison of different flavors, namely those due to 
uncertainties in the unfolding procedure; the relative systematic errors given in the preceding 
section are also applicable, but are common to all three flavor categories; we also list ratios, 
for which these common errors cancel. The flavor unfolding systematic errors were evaluated 
by varying: each diagonal element of the event tagging efficiency matrix en by ±0.01; the 
heavy quark production fractions Rb and R c by the errors on their respective world averages; 
each diagonal bias value B\ by the larger of ±0.005 and ±20% of its difference from unity; 
the photon conversion rate in the simulation by ±15%; and the simulated physics lepton rates 
in light-, c-, and 5-flavor events by ±20%, ±10% and ±5%, respectively. These variations 
correspond to uncertainties on measurements from our data or other experiments except for 
the variation of the bias, which was chosen conservatively to be larger than any change seen 
in a set of comparisons of relevant quantities in data and simulation when selection cuts 
were varied. The unfolding systematic errors are similar in magnitude to (smaller than) the 
statistical errors at low (high) momenta, and are generally dominated by the bias for the 
relevant flavor. There are also substantial contributions from the lepton correction at high 
momentum for and Rg, and from all three bias terms for Rff and R%. 

In figs. 6-8 we show the £ distributions for the three flavor categories, and in fig. 9 we show 
the ratios of production in 6-flavor to light-flavor events and c-flavor to light-flavor events vs. 
x p . At low momentum (high £), there is substantially higher production of charged pions 
in b- and c-flavor events than in light-flavor events, and much higher production of charged 
kaons in 6-flavor events than in light- or c-flavor events; proton production is roughly equal 
at low momentum. As momentum increases, the production of all three charged hadron 
species falls much more rapidly in 6-flavor events than in light-flavor events, and that in 
c-flavor events also drops off sharply at very low £. 

These features are consistent with expectations based on the known properties of e + e~^ 
bb events, namely that a large fraction of the event energy is carried by the leading B and B 
hadrons [33], leaving little energy available to produce nonleading hadrons. The B hadrons 
decay into a large number of lighter particles, including on average 5.5 stable charged hadrons 
[8], which are expected to populate primarily the region 2 < £ < 4. The mixture of particle 
types might be similar to that in light-flavor jets, except that the weak decay chain of the 
B hadron should produce one 'additional' K ± per event, and baryon production might be 
suppressed since a typical baryon-antibaryon pair mass is not small compared with the B 
hadron mass. Similar effects in c-jets result in an intermediate situation: each jet contains 
a charmed hadron with on average about half the beam energy, a lower fraction than for 
5-hadrons, which leaves more energy available for nonleading hadrons than in 6-jets. A D 
meson decay produces an additional kaon that often carries a large fraction of its momentum, 
and fewer additional particles than a B decay. Our results are consistent with those published 
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Pion Production Rates 




Ratios 


<x p > 


uu, dd, ss 


cc 


bb 


cuds 


kuds 


0.0066 


474.0±13.9 


425.5±26.6 


478.1±15.8 


U.oyo±U 


UOo 


1.009±0.019 


0.0088 


467.3±10.5 


440.5±23.2 


488.4±11.9 


n o/i q i n 


U04 


1.045±0.017 


0.0109 


418.2±8.4 


453.8±20.0 


463.7±9.5 


l.UoOztU 


0^^ 
UOO 


1.109±0.018 


0.0131 


375.5±6.9 


409.2±17.2 


432.2±7.7 


i.uyuztu 


UOo 


1.151±0.019 


0.0153 


327.7±5.7 


372.8±14.6 


382.4±6.5 


1 i go i n 
1. looitU 


UOz 


1.167±0.019 


0.0191 


275.8±4.2 


306.4±10.9 


333.3±4.6 


i iii 4_n 


O/I 7 
U4 1 


1.208±0.017 


0.0246 


216.0±3.0 


234.6±8.0 


261.7±3.3 


l.UoOztU 


U44 


1.212±0.018 


0.0301 


171.2±2.2 


197.4±6.2 


214.2±2.7 


1. looztu 


U4o 


1.251±0.017 


0.0356 


140.4±1.9 


155.8±6.0 


175.2±2.3 


i. luyztu 


uou 


1.247±0.016 


0.0411 


116.4±1.5 


132.5±4.3 


145.4±1.9 


i i Qc_i_n 
1. looztU 


U4o 


1.249±0.018 


0.0465 


99.9±1.2 


109.3±3.5 


121.4±1.7 


i.uyoztu 


n/in 

U4U 


1.215±0.018 


0.0521 


85.4±1.0 


92.9±2.9 


103.3±1.5 


i ns8-i-n 
l.UooitU 


Uoo 


1.209±0.019 


0.0576 


72.85±0.89 


77.56±2.48 


89.24±1.36 


l.UOOztU 


0Q7 
Uo I 


1.225±0.020 


0.0630 


64.51±0.79 


68.23±2.17 


75.47±1.21 


l.UOoztU 


UoO 


1.170±0.020 


0.0685 


56.82±0.72 


60.06±1.97 


65.97±1.12 


l.UO ( ztU 


Uoo 


1.161±0.022 


0.0740 


50.84±0.66 


51.76±1.81 


59.39±1.04 


l.UloztU 


Uoo 


1.168±0.023 


0.0795 


45.34±0.61 


45.28±1.67 


52.11±0.97 


u.yyy±u 


OQO 

Uoy 


1.149±0.024 


0.0850 


40.71±0.56 


40.04±1.55 


45.86±0.90 


u.yo4±u 


U41 


1.127±0.025 


0.0931 


34.60±0.40 


33.50±1.12 


38.29±0.65 


u.yooztu 


Uo4 


1.106±0.020 


0.1040 


28.99±0.35 


27.45±0.99 


30.57±0.58 


n Q/l 74-fl 

u.y4 i ztu 


uoo 


1.054±0.021 


0.1150 


24.19±0.31 


22.92±0.87 


24.34±0.51 


n o/i i i o 
u.y4 / ±u 


Uoo 


1.006±0.022 


0.1310 


18.97±0.22 


18.73±0.63 


18.21±0.36 


u.yo i ztu 


uoo 


0.960±0.019 


0.1530 


14.52±0.17 


13.72±0.50 


12.27±0.28 


n Q/i c;-i-n 
u.y40±u 


Uoo 


0.845±0.020 


0.1750 


11.06±0.14 


10.18±0.41 


8.25±0.22 


u.yzu±u 


n/in 

U4U 


0.746±0.021 


0.1970 


8.67±0.12 


7.53±0.34 


5.83±0.18 


U.oOo±U 


U4z 


0.673±0.022 


0.2189 


6.79±0.10 


5.76±0.29 


4.14±0.15 


n s/i Q_i_n 

U.o4oztU 


U40 


0.609±0.023 


0.2410 


5.341±0.085 


4.381±0.235 


2.984±0.128 


n cori-i-n 

U.oZUztU 


U4o 


0.559±0.025 


0.2629 


4.214±0.073 


3.358±0.202 


2.303±0.110 


n 7Q7 1 o 

u. < y t ±u 


UOl 


0.547±0.027 


0.2849 


3.452±0.064 


2.487±0.171 


1.642±0.094 


n 700-1-0 
U. ( ZUztU 


uoz 


0.475±0.028 


0.3068 


2.727±0.056 


1.947±0.148 


1.365±0.085 


O 71 A 4-0 
U. ( 14ztU 


n^7 

UO 1 


0.501±0.032 


0.3338 


2.138±0.042 


1.436±0.108 


0.886±0.063 


O fi79 1 O 
U.O 1 Z±U 


oc;q 

UOo 


0.414±0.030 


0.3666 


1.652±0.036 


0.817±0.087 


0.631±0.052 


0.495±0 


054 


0.382±0.032 


0.3997 


1.164±0.031 


0.614±0.074 


0.490±0.047 


0.527±0 


066 


0.421±0.042 


0.4325 


0.874±0.027 


0.386±0.063 


0.276±0.038 


0.441±0 


074 


0.316±0.045 


U.4001 


0.622±0.024 


0.429±0.061 


0.187±0.033 


0.689±0 


102 


0.301±0.054 


0.5035 


0.441±0.019 


0.206±0.043 


0.111±0.025 


0.467±0 


101 


0.252±0.058 


0.5470 


0.300±0.017 


0.142±0.037 


0.045±0.019 


0.472±0 


128 


0.151±0.062 


0.6083 


0.178±0.010 


0.066±0.021 


0.039±0.010 


0.371±0 


122 


0.222±0.060 


0.7047 


0.081±0.007 


0.003±0.010 


0.011±0.005 


0.043±0 


125 


0.141±0.062 


0.8383 


0.016±0.003 


0.003±0.006 


0.003±0.002 


0.160±0 


356 


0.190±0.150 


Total 


15.294±0.250 


15.783±0.465 


16.841±0.278 





Table 6: Differential production rates (l/N ei , ts )dn 7r ±/dx p of charged pions per Z° decay into 
light (u, d, s), c and b primary flavors. The errors are the sum in quadrature of statistical 
errors and those systematic uncertainties arising from the unfolding procedure; systematic 
errors common to the three flavors are not included in the rates and cancel in the ratios. 
The <x p > values for the three flavor samples are consistent, and have been averaged. 
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Kaon Production Rates 




Ratios 


<iXp> 


mm, dd, ss 


cc 


bb 


cuds 


b.uds 


a ai 
U.Uloo 


27.05±1.27 


30.92±3.86 


30.21±1.99 


1.143±0 


155 


1.117±0.090 


a ai ni 


20.00±0.42 


22.43±1.37 


23.06±0.62 


1.121±0 


078 


1.153±0.040 


n no/is 
U.Uz4o 


19.74±0.40 


22.04±1.27 


22.89±0.60 


1.116±0 


074 


1.159±0.039 


a nom 
U.UoUl 


17.52±0.37 


20.82±1.17 


21.64±0.55 


1.189±0 


076 


1.236±0.041 


a no 
U.UoDO 


16.08±0.37 


16.79±1.15 


21.36±0.53 


1.044±0 


083 


1.328±0.046 


n n/i 1 1 
U.U411 


15.04±0.34 


16.68±1.15 


21.36±0.56 


1.110±0 


086 


1.420±0.048 


U.U40D 


13.54±0.34 


16.46±1.06 


19.90±0.57 


1.215±0 


088 


1.469±0.056 


U.Uozl 


11.87±0.34 


15.81±1.08 


18.91±0.60 


1.332±0 


103 


1.593±0.068 


U.Uo lb 


11.44±0.33 


12.62±0.99 


18.46±0.58 


1.103±0 


094 


1.613±0.069 


U.UboU 


10.64±0.30 


12.24±0.92 


17.43±0.54 


1.151±0 


095 


1.639±0.069 


U.UOOO 


10.24±0.29 


11.42±0.87 


16.92±0.53 


1.115±0 


093 


1.652±0.070 


U.U ( 4U 


9.67±0.29 


10.95±0.85 


15.62±0.52 


1.133±0 


096 


1.616±0.072 


U.U 1 yo 


8.13±0.27 


10.88±0.84 


15.11±0.52 


1.339±0 


114 


1.859±0.090 


U.UooU 


7.98±0.28 


9.62±0.81 


13.18±0.50 


1.206±0 


111 


1.651±0.084 


U.Uyol 


7.00±0.19 


9.84±0.59 


12.43±0.36 


1.406±0 


094 


1.778±0.070 


U. 1U4U 


6.36±0.19 


8.08±0.58 


11.56±0.37 


1.271±0 


100 


1.819±0.080 


n 11 rn 

U.lloU 


5.85±0.20 


8.98±0.63 


9.96±0.38 


1.535±0 


122 


1.704±0.088 


n iQm 
U.iolU 


4.89±0.15 


6.59±0.45 


7.17±0.27 


1.349±0 


104 


1.467±0.072 


U. looU 


3.41±0.17 


5.50±0.51 


4.58±0.29 


1.614±0 


175 


1.343±0.108 


U. 1 1 OU 


2.84±0.22 


5.12±0.68 


4.20±0.36 


1.805±0 


289 


1.480±0.173 


A 1 rWA 

u. iy / u 


2.564±0.082 


3.850±0.245 


2.541±0.126 


1.502±0 


110 


0.991±0.059 


u.zisy 


2.401±0.067 


3.087±0.190 


2.009±0.096 


1.286±0 


090 


0.837±0.047 


U.Z41U 


1.973±0.054 


2.074±0.145 


1.627±0.078 


1.051±0 


081 


0.825±0.046 


u.zozy 


1.643±0.048 


1.960±0.132 


1.116±0.062 


1.193±0 


090 


0.679±0.043 


a oc /i n 
U.Zo4y 


1.481±0.044 


1.681±0.119 


0.830±0.053 


1.135±0 


090 


0.560±0.039 


U.oUoo 


1.211±0.039 


1.368±0.104 


0.640±0.045 


1.129±0 


096 


0.529±0.041 


A QQOQ 

U.oooo 


1.001±0.029 


1.043±0.076 


0.452±0.032 


1.042±0 


084 


0.451±0.034 


U.oooo 


0.746±0.025 


0.874±0.068 


0.337±0.028 


1.171±0 


102 


0.451±0.040 


a Qno7 

u.oyy l 


0.666±0.023 


0.600±0.058 


0.245±0.024 


0.900±0 


094 


0.367±0.038 


U.4oZo 


a ct^nxn aoo 

u.ooy±u.uzz 


a /incxn a £A 
U.4Uo±U.UoU 


A 1 A A_1_A AOA 

u.i4y±u.uzu 


0.730±0 


094 


U.ZOO±U.UoO 


0.4651 


0.426±0.020 


0.408±0.050 


0.108±0.018 


0.957±0 


126 


0.253±0.043 


0.5035 


0.363±0.016 


0.243±0.037 


0.057±0.012 


0.669±0 


106 


0.156±0.035 


D ^zL7D 

U.O^J: 1 U 


0.261±0.015 


0.173±0.034 


0.061±0.013 


0.663±0 


136 


0.233±0.050 


0.6083 


0.183±0.010 


0.064±0.020 


0.012±0.005 


0.351±0 


109 


0.066±0.030 


0.7047 


0.079±0.007 


0.009±0.011 


0.002±0.003 


0.113±0 


141 


0.026±0.035 


0.8383 


0.008±0.004 


0.008±0.008 


-.001±0.001 


0.989±1 


126 


-.098±0.102 


Total 


1.869±0.062 


2.273±0.093 


2.377±0.080 





Table 7: Differential production rates (1/N e „ ts )dnx± / dx p of K ± mesons per Z° decay into 
light, c and b primary flavors. 
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Proton Production Rates 






Ratios 


<x p > 


— 7 7 — 

mm, act, ss 


cc 


7 7 

bb 




cuds 


7 7 

b:uds 


U.UiDo 


13.98±0.99 


13.28±2.94 


13.79±1.49 


U 


950±0 


226 


0.987±0.127 


n ni 01 


17.63±0.58 


15.22±1.89 


17.93±0.78 


u 


863±0 


120 


1.017±0.056 


U.UZ40 


13.42±0.60 


13.32±1.86 


16.41±0.93 


n 
U 


992±0 


152 


1.223±0.089 


n n^m 
U.UoUl 


10.57±0.36 


9.60±1.16 


12.11±0.52 


n 
u 


909±0 


122 


1.146±0.063 


U.Uooo 


9.98±0.31 


11.64±1.01 


10.32±0.42 


1 


166±0 


118 


1.034±0.052 


n n/ii 1 

U.U41 1 


8.37±0.26 


10.07±0.87 


9.52±0.40 


1 


203±0 


118 


1.138±0.059 


U.U40O 


7.33±0.24 


8.10±0.76 


7.72±0.36 


1 


105±0 


117 


1.054±0.059 


n n^oi 


7.79±0.23 


6.09±0.72 


6.86±0.32 


n 
u 


781±0 


099 


0.879±0.049 


U.UO ( 


6.62±0.22 


6.54±0.67 


6.19±0.32 


n 
u 


988±0 


111 


0.935±0.057 


U.UOoU 


5.88±0.20 


6.36±0.63 


4.96±0.29 


1 


082±0 


117 


0.845±0.056 


U.UOoO 


5.39±0.19 


4.62±0.59 


4.82±0.29 


n 
u 


857±0 


115 


0.895±0.062 


n n7/in 
U.U 1 4U 


5.22±0.19 


4.43±0.58 


4.57±0.29 


n 
u 


848±0 


117 


0.875±0.064 


u.u / yo 


4.42±0.18 


4.08±0.55 


4.07±0.29 


n 
u 


924±0 


133 


0.920±0.076 


n ns c;n 
U.UoOU 


4.44±0.19 


3.67±0.55 


3.82±0.29 


u 


827±0 


131 


0.860±0.074 


U.Uyol 


3.65±0.13 


4.07±0.42 


3.29±0.21 


1 
1 


116±0 


123 


0.901±0.066 


n 1 n/in 

U. 1U4U 


3.11±0.13 


2.98±0.41 


2.68±0.22 


u 


960±0 


139 


0.861±0.079 


n 1 1 en 
U. IIOU 


2.73±0.15 


2.30±0.43 


2.24±0.23 


n 
u 


844±0 


166 


0.819±0.096 


n 1 qi n 
U. lolU 


2.15±0.12 


2.39±0.36 


1.84±0.19 


1 


110±0 


180 


0.856±0.100 


U. 100U 


1.83±0.16 


1.72±0.50 


1.91±0.27 


n 
U 


940±0 


287 


1.048±0.174 


U. 1 ( OU 


1.84±0.24 


0.31±0.71 


1.25±0.36 


n 
U 


171±0 


392 


0.683±0.219 


n 1 o7n 
u. iy < u 


0.905±0.078 


0.561±0.235 


0.867±0.121 


U 


619±0 


265 


0.958±0.157 


n 01 so 


1.065±0.054 


0.978±0.163 


0.739±0.072 


u 


918±0 


161 


0.695±0.077 


n o/i 1 n 
U.241U 


0.822±0.044 


0.907±0.136 


0.645±0.060 


1 
1 


104±0 


176 


0.784±0.085 


u.zozy 


0.762±0.038 


0.652±0.116 


0.392±0.047 


u 


855±0 


159 


0.515±0.067 


u.zo4y 


0.628±0.033 


0.572±0.101 


0.252±0.038 


n 
U 


911±0 


169 


0.401±0.064 


U.0UO0 


0.486±0.029 


0.494±0.089 


0.266±0.035 


1 


016±0 


195 


0.547±0.079 


U.0000 


0.446±0.022 


0.454±0.069 


0.146±0.022 


1 
1 


016±0 


165 


0.327±0.052 


u.0000 


0.306±0.018 


0.314±0.054 


0.102±0.018 


1 
l 


026±0 


190 


0.333±0.061 


n QQQ7 

u.oyy < 


0.230±0.015 


0.170±0.041 


0.020±0.012 


u 


741 ±0 


186 


0.085±0.053 


U.4OZ0 


n 1 Q7_i_n ni q 
u. iy t ztu.uio 


U. lUoztU.Uoo 


U.Uo4±U.Ull 


n 
u 


522±0 


171 


n 1 71 -i~n n^c; 

U. 1 ( IztU.UOO 


0.4651 


0.145±0.011 


0.064±0.028 


-.004±0.007 





445±0 


196 


-.026±0.048 


0.5035 


0.108±0.008 


0.015±0.018 


0.016±0.006 





142±0 


166 


0.150±0.058 




0.070±0.006 


0.044±0.017 


-.003±0.003 





626±0 


243 


-.043±0.049 


0.6083 


0.036±0.003 


0.007±0.007 


0.004±0.002 





186±0 


193 


0.107±0.061 


0.7047 


0.013±0.002 


0.015±0.005 


-.001±0.001 


1 


150±0 


434 


-.089±0.045 


0.8383 


0.003±0.001 


0.001±0.002 


0.000±0.000 





221±0 


620 


0.159±0.160 


Total 


1.008±0.038 


0.930±0.056 


0.909±0.037 





Table 8: Differential production rates (l/N et , ts )dn p /p/da;p of p and p per Z° decay into light, 
c and b primary flavors. 
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8-2003 ?— InCx ) 

8675A5 m W 

Figure 6: Distributions of £ in light-flavor events for (a) pions (circles) and (b) kaons (squares) 
and protons (diamonds), compared with the predictions of the three models. 
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Figure 7: Distributions of £ in c-flavor events for (a) pions (circles) and (b) kaons (squares) 
and protons (diamonds), compared with the predictions of the three models. 
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Figure 8: Distributions of ^ in 6-flavor events for (a) pions (circles) and (b) kaons (squares) 
and protons (diamonds), compared with the predictions of the three models. 
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8675A8 X P= 2 P/ECM 

Figure 9: Ratios of pion (top), kaon (middle) and proton (bottom) production in 6-flavor 
events to that in light-flavor events (left) and in c-flavor: light-flavor events (right). Some 
bins have been combined for clarity. Also shown are the predictions of the three models. 
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<x p > 



0.0066 
0.0088 
0.0109 
0.0131 
0.0153 
0.0191 
0.0246 
0.0301 
0.0356 
0.0411 
0.0465 
0.0521 
0.0576 
0.0630 
0.0685 
0.0740 
0.0795 
0.0850 
0.0931 
0.1040 
0.1150 
0.1310 
0.1530 
0.1750 
0.1970 
0.2189 
0.2410 
0.2629 
0.2849 
0.3068 
0.3338 
0.3666 
0.3997 
0.4325 
0.4651 
0.5035 
0.5470 
0.6083 
0.7047 
0.8383 



Charged Particle Production Rates 



uu, dd, ss 



7.5 
6.3 
5.8 
4.2 
3.2 
2.3 
2.0 
1.6 
1.3 
1.1 
1.0 



507.8± 11.9 
505.2± 9.2 
465.3± 
421. 9± 
371. 7± 
315.5± 
250.5± 
200.3± 
167.3± 
140.4± 
121. 2± 
105.5± 
91. 2± 
81.29± 0.89 
72.69± 0.81 
65.92± 0.76 
58.06± 0.70 
53.26± 0.66 
45.37± 0.48 
38.55± 0.43 
32.84± 0.41 
26.05± 0.28 
19.79± 0.29 
15.75± 0.35 
12.16± 0.17 
10.27± 0.14 
8.14± 0.11 
6.62± 0.10 
5.565±0.087 
4.428±0.076 
3.588±0.057 
2.706±0.049 
2.062±0.042 
1.631±0.037 
1.193±0.034 
0.912±0.026 
0.632±0.023 
0.398±0.015 
0.172±0.011 
0.027±0.005 



cc 



9.5 
7.3 
6.9 
5.2 
4.2 
3.6 
3.0 



468.7± 25.3 
485.4± 22.6 
507.2± 20.2 
464.0± 17.6 
422.9± 16.9 
349.1± 12.7 
274. 1± 
231. 5± 
187.5± 
162.3± 
136.6± 
117.3± 
99. 1± 
89.21± 2.65 
78.25± 2.42 
69.26± 2.23 
62.25± 2.06 
55.28± 1.93 
49.14± 1.42 
40.11± 1.29 
35.64± 1.22 
28.94± 0.83 
21.99± 0.86 
16.51± 1.03 
12.69± 0.50 
10.41± 0.40 
7.86± 0.32 
6.37± 0.28 
5.060±0.241 
4.080±0.210 
3.123±0.156 
2.141±0.126 
1.472±0.106 
0.952±0.090 
0.935±0.085 
0.485±0.061 
0.372±0.055 
0.140±0.031 
0.027±0.016 
0.011±0.010 



bb 



8.0 
6.3 
6.7 
4.4 
3.4 
2.6 
2.1 
1.9 
1.7 
1.6 
1.4 
1.3 



546. 1± 14.2 
558.5± 10.3 
531.9± 
490.8± 
436.5± 
382.8± 
308.6± 
254.8± 
213.2± 
182. 1± 
154.5± 
134.3± 
118.6± 
102.4± 
91.92± 1.22 
83.63± 1.14 
75.06± 1.10 
66.58± 1.04 
57.31± 0.74 
47.80± 0.71 
39.19± 0.67 
29.54± 0.48 
20.69± 0.48 
15.36± 0.55 
10.65± 0.26 
8.06± 0.20 
6.28± 0.16 
4.69± 0.13 
3.490±0.110 
2.935±0.098 
2.041±0.068 
1.534±0.058 
1.111±0.051 
0.736±0.042 
0.510±0.037 
0.330±0.027 
0.188±0.022 
0.089±0.012 
0.017±0.006 
0.003±0.002 



Ratios 



cuds 



0.923±0.053 
0.961±0.050 
1.090±0.052 
1.100±0.050 
1.137±0.054 
1.107±0.049 
1.094±0.046 
1.156±0.044 
1.121±0.050 
1.156±0.044 
1.127±0.041 
1.112±0.040 
1.087±0.038 
1.097±0.037 
1.077±0.037 
1.051±0.038 
1.072±0.040 
1.038±0.040 
1.083±0.035 
1.040±0.037 
1.085±0.041 
1.111±0.035 
1.111±0.048 
1.048±0.071 
1.044±0.046 
1.013±0.044 
0.965±0.044 
0.961±0.047 
0.909±0.048 
0.922±0.052 
0.870±0.048 
0.791±0.051 
0.714±0.055 
0.584±0.058 
0.783±0.076 
0.532±0.070 
0.588±0.091 
0.353±0.081 
0.158±0.096 
0.412±0.386 



kuds 



1.075±0.014 
1.105±0.013 
1.143±0.014 
1.163±0.014 
1.174±0.020 
1.213±0.017 
1.232±0.018 
1.272±0.017 
1.275±0.016 
1.297±0.017 
1.275±0.018 
1.273±0.019 
1.301±0.020 
1.260±0.020 
1.265±0.021 
1.269±0.022 
1.293±0.024 
1.250±0.024 
1.263±0.020 
1.240±0.022 
1.194±0.025 
1.134±0.022 
1.045±0.029 
0.975±0.041 
0.876±0.024 
0.785±0.022 
0.772±0.022 
0.709±0.023 
0.627±0.022 
0.663±0.025 
0.569±0.021 
0.567±0.024 
0.539±0.027 
0.451±0.028 
0.427±0.033 
0.362±0.032 
0.298±0.037 
0.224±0.031 
0.099±0.034 
0.103±0.095 



Table 9: Differential production rates (l/~N evts )dn chg /dx p of stable charged particles per Z° 
decay into light, c and b primary flavors. 
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previously [28] and considerably more precise. 

Also shown in figs. 6-9 are the predictions of the three models; they all reproduce the 
observed flavor dependence qualitatively. In the case of the light flavors, problems with the 
models very similar to those seen in the preceding section for the flavor inclusive sample 
are observed, indicating problems in the modelling of hadronization, and not just in that 
of heavy hadron decay. Some of the discrepancies are larger in the light flavors than in all 
flavors; in particular the structure at small £ in the HERWIG model is now very pronounced 
for all three particle types. 

In c-flavor events, all models predict the pion yield within a few percent, but the spectra 
are slightly too hard. JETSET again predicts the shape of the proton spectrum reasonably 
well, but is slightly high on the amplitude; UCLA and HERWIG describe the data at large £, 
but fall well below the data at small £. JETSET and UCLA predict a kaon spectrum in which 
those from leading D hadron decays peak at £ ps 2, and those from the hadronization of the 
remaining jet form a broad shoulder at larger £ values. They are qualitatively consistent with 
the data, though the data prefer a lower shoulder. HERWIG predicts a quite different shape 
that is inconsistent with the data, probably reflecting known problems in the modelling of 
heavy hadron production and decay. 

In 6-flavor events, JETSET and UCLA describe the pion and kaon spectra reasonably 
well, predicting a few too many pions in the range 2.5 < £ < 4, and too many kaons for 
£ > 3.5; both also predict the shape of the proton spectrum well, but UCLA is too low in 
amplitude and JETSET is consistent with the data. Again, HERWIG predicts very different 
spectra, all of which are inconsistent with the data at large £ values. 

The rate for all charged tracks in each of the flavor samples was derived by summing 
the rates for the three charged hadron species and adding the simulated contribution from 
leptons. These rates and their ratios are listed in table 9; no plots are shown since the features 
correspond to those of the pion data and models in figs. 6-9. Our data are consistent with 
those published by DELPHI [28] and OPAL [10], and substantially more precise. 

7 Comparison with QCD Predictions 

We tested the predictions of QCD in the Modified Leading Logarithm Approximation (MLLA), 
combined with the ansatz of Local Parton-Hadron Duality (LPHD), by fitting Gaussian and 
distorted Gaussian functions to our measured £ distributions. Examples of such fits are 
shown in figs. 10 and 11. In each case, we first fitted a simple Gaussian to each set of points 
within a region of £ at least 0.5 units in size centered on the approximate peak position. 
We chose the largest such symmetric range for which the confidence level of the x 2 of the 
fit exceeded 0.01, and then added points on one side (in all cases at lower £ values) as long 
as the confidence level remained above 0.01. We found that the Gaussian function could 
describe the data over a symmetric range of at least 0.8-1.6 units about the peak position, 
consistent with one prediction of MLLA QCD; the fit could be extended toward lower £ by 
as much as 1.7 additional units, although it must be noted that our errors on the kaon and 
proton spectra are rather large in this range. The largest fittable ranges are given in table 10 
and the corresponding fitted functions are shown on figs. 10 and 11. 

We next introduced a skewness term into the function (G + ) and repeated the above 
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procedure. We found that the symmetric range could be extended in only some cases, and 
by at most 0.6 units (see table 10). In all-flavor and light-flavor events there was always an 
increase in the fittable range on either the low or high £ end, and the fitted skewness values 
were small. In c- and 6-flavor events, however, there were some cases in which no increase in 
the fittable range could be obtained, and others in which the range could be increased but 
the skewness value increased rapidly. 

The addition of a kurtosis term (G ++ ) had similar results (see table 10); the fittable 
range was increased greatly, in many cases to the entire measured range, however both the 
skewness and kurtosis values became large for c- and 6-flavor events. The resulting G ++ 
functions are shown on figs. 10-11 over the entire £ range; the large distortions are evident 
in the c- and 6-flavor events for kaons and protons. The MLLA prediction that a Gaussian 
with small distortion terms should describe the data over a range substantially larger than 
one unit about the peak position holds for the light flavors, but does not hold for heavy 
flavors; this might be expected since the calculation assumes massless partons. 

The peak £* of the £ distribution is predicted to decrease exponentially with increasing 
particle mass. Following convention, we took the mean of the fitted Gaussian over a range of 
one unit about the peak as an estimate of £*, and, in addition to statistical and experimental 
systematic errors, we considered a variation of the fit range. A fit was performed to each 
set of contiguous points with a smallest (largest) £ value between 0.75 and 1.25 units below 
(above) the peak position. Half the difference between the highest and lowest of the fitted 
means was taken as an estimate of the systematic uncertainty due to the fit range. The 
resulting values of are listed in table 10, where a considerable flavor dependence is seen. 

The values in all-flavor and light-flavor events are shown in fig. 12 along with previous 
results [7, 12, 27, 28]. As observed previously, the £* values for the measured hadron species 
in all-flavor events do not fall on a single trajectory as a function of mass, although roughly 
parallel trajectories can be postulated for mesons and baryons. For light flavors, the £* value 
for pions is quite similar to that for all flavors, and those for the other mesons (baryons) tend 
to be higher (lower), as would be expected if the MLLA prediction is correct for primary 
hadrons, but the observed distributions are distorted by secondary particles. However it 
is still not possible to draw a single trajectory through all the light-flavor points and this 
apparent deficiency in the the theory remains; it is possible that the exclusion of additional 
non-primary hadrons might remove it. 

Our previous results [12] were used in conjunction with data from other experiments to 
test the perturbative QCD calculations of the E C m dependence of the rates at high-x p [5]. 
The results presented here can be used to make a slightly more precise test; however more 
precise data at lower and/or higher E CM are needed to confront these calculations in a more 
meaningful way. 

8 Total Production Rates 

We have integrated our differential production rates over their respective measurement 
ranges, taking into account the bin-to-bin correlations in the systematic errors. These inte- 
grated rates per event are listed in tables 2-5 and 6-8; the errors are dominated by overall 
normalization uncertainties corresponding to the uncertainty in our track reconstruction ef- 
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G 


[ax. fit rari£ 
G + 


C++ 


Peak position 
£*±stat. ±syst. ±fit 


All 

tt ± : uds 

c 

b 


1.62-4.74 
1.40-4.85 
1.30-4.60 
1.81-5.21 


1.88-5.21 
1.55-5.21 
1.47-5.21 
1.81-5.21 


0.92-5.21 
0.64-5.21 
0.87-5.21 
0.00-5.21 


3.752±0.004±0.008±0.010 
3.787±0.008±0.027±0.014 
3.726±0.014±0.090±0.039 
3.684±0.008±0.005±0.001 


All 

K ± : uds 
c 
b 


0.92-4.18 
0.55-4.25 
0.87-3.61 
0.72-3.61 


0.49-4.18 
0.24-4.25 
0.00-4.25 
0.00-4.25 


0.34-4.18 
0.24-4.25 
0.00-4.25 
0.00-4.25 


2.549±0.018±0.042±0.062 
2.592±0.032±0.091±0.021 
2.412±0.030±0.018±0.012 
2.587±0.014±0.017±0.004 


All 

pp: uds 
c 
b 


0.84-4.18 
0.64-4.25 
0.72-4.25 
0.97-4.25 


0.76-4.18 
0.42-4.25 
0.72-4.25 
0.97-4.25 


0.14-4.18 
0.00-4.25 
0.00-4.25 
0.00-4.25 


3.084±0.101±0.056±0.031 
2.859±0.076±0.028±0.026 
3.079±0.113±0.023±0.068 
3.513±0.100±0.254±0.011 



Table 10: The widest ranges in £ (see text) over which a Gaussian function (G) was able to 
describe the data alone, and with the addition of skewness (G + ) and kurtosis terms (G ++ ). 
Peak positions £* from the Gaussian fits described in the text; the errors are statistical, 
experimental systematic, and due to variation of the fit range. 
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Figure 10: Measured £ distributions for all charged particles (top) and pions (bottom) in all-, 
light-, c- and 6-flavor events. Some bins have been combined for clarity. The solid (dashed) 
lines represent the results of the maximal Gaussian (G) and G ++ fits described in the text. 
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Figure 11: Measured £ distributions for charged kaons (top) and protons (bottom) in all-, 
light-, c- and 6-flavor events. Some bins have been combined for clarity. The solid (dashed) 
lines represent the results of the maximal Gaussian (G) and G ++ fits described in the text. 
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Figure 12: Peak positions as a function of hadron mass in (a) all- and (b) light-flavor 
events, along with previous results from experiments at the Z°. 



ficiency. In order to quote total rates, we must extrapolate into the unmeasured regions, 
£ > 5.21 for inclusive charged particles and pions, and £ > 4.25 for kaons and protons. As 
can be seen in figs. 5-8, and 10-11, this is a non-trivial effect, and none of the models or 
functions provides an ideal estimate of the accepted fraction. However, a set of four esti- 
mates derived from the three models and the G ++ function might be considered to cover 
a reasonable range of possibilities. For the c- and 6-flavor events, we used in addition the 
generator used for our detector simulation [26], a version of JETSET tuned to the world's 
data on D and B hadron production and decay. 

From the pions and protons generated using each of the models, we calculated the fraction 
that were generated in the range of our measurement, as the predictions, when normalized 
to the data, provided reasonable descriptions of the shape. For the charged kaons (figs. 5- 
8), all models predict spectra that are too soft; we therefore scaled the predictions along 
the £ axis until the best agreement with the data was found. This procedure changed the 
acceptance correction by a few percent for b- and c-flavor events, and by as much as 12% for 
the HERWIG model in light-flavor events. The fitted G ++ function was integrated to obtain 
another estimate of the accepted fraction in those cases where it gave an acceptable \ 2 over 
the entire measured range. Otherwise, it was used to calculate a fraction above the lower 
bound given in table 10, and the data below that bound were added to obtain an overall 
fraction. 

For each hadron species and flavor category these four or five estimates of the fraction 
were found to be similar, with a typical rms value of about 1% relative. We took their average 
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Coverage 


Yield/Event 


Difference 


All 

n ± : uds 
c 
b 


0.925±0.007 
0.923±0.008 
0.931±0.013 
0.929±0.007 


17.007±0.209 
16.579±0.304 
16.954±0.556 
18.136±0.330 


c—uds 
b—uds 
b-c 


0.375±0.587 
1.557±0.232 
1.182±0.498 


All 

K ± : uds 
c 
b 


0.941±0.010 
0.934±0.007 
0.937±0.005 
0.947±0.006 


2.203±0.071 
2.000±0.068 
2.427±0.100 
2.510±0.086 


c—uds 
b—uds 
b-c 


0.427±0.074 
0.510±0.037 
0.083±0.075 


All 

pp: uds 
c 
b 


0.933±0.003 
0.921±0.010 
0.899±0.038 
0.906±0.020 


1.054±0.035 
1.094±0.043 
1.034±0.077 
1.004±0.046 


c—uds 
b—uds 
b-c 


-0.060±0.074 
-0.091±0.034 
-0.031±0.074 


uds 

All Chg.: c 
b 


0.927±0.007 
0.940±0.015 
0.933±0.010 


20.048±0.316 
21.096±0.653 
23.098±0.378 


c—uds 
b—uds 
b-c 


1.048±0.718 
3.050±0.311 
2.002±0.643 



Table 11: Estimated fractions of the total production rates contained within the range of 
the measurements; corrected yields of charged hadrons per Z° decay into events of each 
flavor category. The rightmost columns show differences between flavors, for which some 
uncertainties cancel. 
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as our central value of the fraction in each case, and took their rms as the uncertainty due to 
the extrapolation procedure. These values along with the corrected total rates are listed in 
table 11; also shown are differences between the three flavor categories, for which some of the 
uncertainties cancel. The results were found to be consistent with previous measurements; 
the precision is similar for all-flavor events, and substantially improved for light- and 6-flavors; 
ours are the only measurements for c-flavor events. 

Roughly 10% more pions are produced in 6-flavor events than in c- or light-flavor events; 
roughly 20% more kaons are produced in both b- and c-flavor events than in light-flavor 
events; roughly 10% fewer protons are produced in b- than in light-flavor events. The total 
charged multiplicities and differences between flavors are consistent with previous, dedicated 
measurements [9, 10]; they have comparable precision with different systematic error sources. 



9 Leading Particle Effects 



We extended these studies to look for differences between particle and antiparticle production 
in light quark (as opposed to antiquark) jets, in order to address the question of whether 
e.g. a primary -u-initiated jet contains more hadrons that contain a valence -u-quark (e.g. 
7r + , K + , p) than hadrons that do not (e.g. n~ , K~ , p). To this end we used the light quark- 
and antiquark-tagged hemispheres described in section 3. 

We measured the differential production rates per light quark jet 



d 



m = 1 !/ 



N(q -> h) + N(q -> h) 
N(q ^h)+ N(q -> h) 



(2) 
(3) 



where: q and q represent light-flavor quark and antiquark jets respectively; N evts is the 
total number of events in the sample; h represents n + , K + or p, and h the corresponding 
antihadron. Then, for example, N(q — > h) is the number of hadrons of species h in light 
quark jets. This formulation assumes CP symmetry, i.e. N(q — > h) — N(q — > h), which was 
found to be satisfied in the data in all cases. 

The charged hadron analysis was repeated on the sample of positively charged tracks in 
the quark-tagged jets and negatively charged tracks in the antiquark-tagged jets, yielding 
measured values of R^+, R q K +, and in the tagged samples. The same procedure applied 
to the remaining tracks yielded R q K ~, and Rp. 

The decays of the leading heavy hadrons in simulated heavy flavor background events give 
rise to differences between hadron and antihadron production in the quark-tagged sample 
over the entire x p range, including an effect in high-momentum pions of opposite sign to 
that seen in the light-flavor data (below). It is essential to understand and/or suppress 
this contribution; our simulation has been tuned to the available data on charmed and 
bottom hadron decays, and our light-flavor tag reduces this background to 9% cc and 2% 
bb events. At this level the simulated contribution to any difference is well below the other 
uncertainties. The simulated contribution to each rate was applied as a correction, yielding 
differential production rates per light-quark-tagged jet. For each hadron species, differential 
production rates in light quark jets were then extracted by correcting for the light-tag bias 
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Table 12: Differential production rates R\ = (l/2N evts )drih/ dx p for positively and negatively 
charged pions h = n + , tt~ in light (u, d and s) quark jets from hadronic Z° decays, along 
with the normalized difference D n - = (R^- — -^+)/(-^- + ^+)- The errors are the sum 
in quadrature of statistical errors and those systematic errors arising from the light quark 
tagging and unfolding procedure. 



(see sec. 6) and unfolding for the effective quark (vs. antiquark) purity. This purity depends 
on polar angle, and hence on the acceptance of the CRID as well as the event selection 
discussed in section 3; its average value was estimated from the simulation to be 0.72 for the 
selected track sample. 

The measured differential production rates per light quark jet are listed in tables 12-13 
and shown in fig. 13; as for the flavor dependent results (sec. 6), the errors given are the sum 
in quadrature of the statistical error and those systematic errors arising from the tagging 
and correction procedures. The latter include variation of: the event tagging efficiencies 
and biases as described in section 6; the electroweak parameters Rb, R c , Ab and A c by the 
errors on their respective world average values; the effective quark purity by ±0.01; the 
sum of h and h rates in c- and 6-flavor events by a smooth parametrization of the errors in 
tables 6-8; and their difference by ±20% of itself to cover the uncertainty in the electron 
beam polarization and statistical error on the simulated purities. The systematic errors are 
small compared with the statistical errors, and are typically dominated by the uncertainty 
on the effective quark purity. 
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Table 13: Differential production rates R q h = (l/2N evts )dnh/ dx p for positively and negatively 
charged kaons and protons in light (w, d and s) quark jets from hadronic Z° decays, along 
with the normalized differences Dh = {R^ — Rj^)/{R q h + The errors are the sum in 

quadrature of statistical errors and those systematic errors arising from the light quark 
tagging and unfolding procedure. 
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Figure 13: Fully corrected production rates for positively and negatively charged hadrons 
in light quark jets. The error bars include statistics and those systematic uncertainties that 
affect the difference between the charges. 
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In all cases the hadron and antihadron rates are consistent at low x p . For kaons (protons), 
significant differences appear for x p above about 0.15 (0.2) units and grow with increasing 
x p . For pions, differences are smaller, but significant for x p above about 0.3, also appearing 
to grow with x p . It is convenient to show these results in the form of the difference between 
hadron h and antihadron h production normalized by the sum: 

The common systematic errors cancel explicitly in this variable, which is shown for the 
hadrons h = tc~,K~, p in fig. 14 and listed in tables 12-13. A value of zero corresponds 
to equal production of hadron and antihadron, whereas a value of +(-)l corresponds to 
complete dominance of (anti)particle production in light quark jets. 

The results for the protons afford the most straightforward interpretation. Since baryons 
contain valence quarks and not antiquarks, the positive values of D p for x p > 0.2 are clear 
evidence for the production of leading protons. The data are consistent with a steady 
increase with x p to a plateau of D p rs 0.5 for x p > 0.5, although the errors on the highest 
x p points are quite large. For x p < 0.1 the data are consistent, within common systematic 
errors, with equal production of baryons and antibaryons, however it must be noted that 
the contribution from non-leading hadrons is very high in this region and we cannot exclude 
that some leading baryons are produced at low x p . 

The interpretation of our results for mesons is more complicated, since a meson contains 
one valence quark along with one valence antiquark. In SM Z° decays all primary down- 
type quarks are produced equally and with the same forward-backward asymmetry, so that 
if a leading neutral hadron such as K*° (sd) were produced equally in s and d jets then 
we would measure D^*o to be zero. We previously measured [12] a significantly positive 
value (see fig. 14) at high x P) indicating both that there is leading K*° production and that 
more leading K*° are produced in s jets than in d jets. This is an expected consequence of 
strangeness suppression in the hadronization process. That is, it is expected to be less likely 
for an ss pair to pop out of the vacuum and the s to pair up with an initial d than it is for 
a dd pair to pop out and the d to pair up with the initial s. 

In the case of charged mesons such as n~ (du), there is a nonzero dilution of leading parti- 
cle effects because of the different Z° branching ratios and forward-backward asymmetries of 
up- and down-type quarks. Assuming SM couplings to the Z°, equal production of leading 7r + 
in -u-jets and 7r~ in d-jets, and no contribution from s-jets, we calculate a dilution factor for 
this analysis of 0.27. That is, we would expect to observe Dfl = 0.27 = 0.27 D^+. 

Our measured D n - are significantly positive at high x p , indicating leading pion produc- 
tion. The data give no information on the relative contributions of the three light flavors, 
but are consistent with 0.27-D p at all x p , and hence with the notion that leading pion effects 
in u and d jets are of similar strength to those of leading baryons in all light quark jets. Our 
measured D K - are consistently positive for x p > 0.15, and well above both 0.27_D p and D n - 
at high x p . This indicates both production of leading charged kaons and greater production 
of leading K~ in s-jets than in -u-jets. 

The predictions of the three models are also shown in fig. 14. For protons, the HERWIG 
prediction drops below zero for 0.05 < x p < 0.25, then rises rapidly to unity at higher x p ; 
this structure is inconsistent with our data. The UCLA prediction also rises to unity at high 
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Figure 14: Normalized differences between hadron and antihadron production in light quark 
jets, compared with the predictions of the three models. Some bins have been combined 
for clarity. Our previous results for A (K*°) are overlaid on the proton (kaon) plot; the 
corresponding model predictions are similar. 
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Xp, becoming inconsistent with the data for x p >0.25. The JETSET prediction is consistent 
with the data. For the mesons, all models predict positive values of D n - and D K - even 
at very low x p ; they are all consistent with our data in that region, but we are not able to 
resolve the difference from zero. At high x p , all three predictions for D n - lie roughly 2er 
below the data, and the measured D K - favor the JETSET prediction slightly over those of 
HERWIG and UCLA. 

10 Summary and Conclusions 

We have improved our measurements of the production of the charged hadron species n ± , 
K ± and p/p, as well as that of inclusive stable charged particles, in hadronic Z° decays, tak- 
ing advantage of the sample of 400,000 decays recorded with the upgraded vertex detector 
to reduce both the statistical and systematic errors substantially. The SLD Cherenkov Ring 
Imaging Detector enabled the clean and efficient identification of charged tracks over a wide 
momentum range, yielding precise measurements of their total and differential production 
rates. Our measurements in flavor-inclusive hadronic Z° decays are consistent with, com- 
plementary to, and in some regions more precise than previous measurements. Deficiencies 
in popular hadronization models have been confirmed. 

The precision of the vertex detector allowed us to isolate very high-purity light-, c- and 
6-tagged event samples, and to compare production characteristics of the hadron species in 
light-, c- and 6-flavor events. Significant differences between flavors were found, consistent 
with expectations based on the known properties of B and D hadron production and decay. 
The known problems with hadronization models were all present in the light-flavor events, 
confirming that they are indeed in the hadronization itself and not just artifacts of heavy 
hadron modelling. Additional problems with the models were observed in c- and 6-flavor 
events. 

The shape of the £ = — In x p distribution near its peak for each hadron species in events 
of each flavor is consistent with the Gaussian form predicted by MLLA QCD+LPHD. For c- 
and 6-flavor events, however, the Gaussian form cannot accommodate the data over a wider 
range without the addition of very large distortion terms. In light-flavor events, the Gaussian 
with small additional distortion terms is able to describe the data over a wider range. The 
peak positions £* for each hadron species in light-flavor events are more consistent with a 
monotonic dependence on hadron mass than those in flavor-inclusive events. Our data are 
thus consistent with the predictions of MLLA QCD for the light flavors, and indicate that 
the presence of heavy hadrons distorts the observed spectra. 

Using the large forward-backward asymmetry induced by the polarized SLC electron 
beam to separate light quark from light antiquark hemispheres, we have compared hadron 
and antihadron production in light quark jets. A large excess of proton over antiproton 
production at high x p is direct evidence for the production of leading baryons, i.e. baryons 
that carry the quantum numbers of the quark that initiated the jet. A similarly large 
excess of K~ over K + production indicates not only leading kaon production but also that 
leading kaons are produced more often from initial s quarks than from initial u quarks. 
A smaller excess of ir~ over 7r + production was also observed. This is expected if leading 
pions account for a large fraction of high momentum pion production but with equal 7r~ 
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rates from initial d and u jets: uu and dd events then contribute to the observed difference 
with opposite signs, but since down-type quarks are produced more often and with a higher 
electroweak asymmetry than up-type quarks, a net positive difference is observed. These 
data provide unique and stringent tests of hadronization models. All models tested were 
able to reproduce the pion and kaon data, though the latter favor the JETSET model over 
the other two; JETSET is also consistent with the proton data, but the other two models 
predict values much higher than the data for x p >0.3. 
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